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(N ■ Abstract In this paper, we demonstrate that interleaved sampling techniques can 

. be used to characterize the Hamiltonian of a qubit and its environmental decoher- 

ence rate. The technique offers a significant advantage in terms of the number of 
measurements that are required to characterize a qubit. When compared to the 
standard Nyquist-Shannon sampling rate, the saving in the total measurement 
time for the interleaved method is approximately proportional to the ratio of the 
sample rates. 
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The qubit is the principal building block for most quantum information processing 
schemes. The ability to manipulate the quantum state of one or more qubits to 
(N ■ perform useful quantum computing operations is reliant on a number of factors. 

In this paper we consider two: the estimation of the qubit Hamiltonian and the 
characterization of environmental decoherence via a sequence of projective mea- 
surements. The use of projective measurements to characterize a qubit is not new, 
it has been demonstrated by a number of groups for single qubits and two qubits 
^ (e.g. [Il[2l[3]) and methods for interpreting the measurements to reconstruct the 

^— V . Hamiltonian and the environmental decoherence rate have been considered (e.g. [H 

o : ElEl)- 

In this paper, we propose a characterization method based on classical sam- 
pling methods which could be a more efficient way to determine the dynamical 
parameters of an experimental qubit. Specifically, we demonstrate that, when the 
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decoherence in the qubit is low enough for the frequency components of the co- 
herent osciUations to be confined to a relatively small band of frequency space, 
signal reconstruction techniques based on two sets of interleaved measurements 
provide robust estimates of the relevant parameters with significantly fewer mea- 
surements. The reduction in the number of measurements is very important in 
this case because the use of projective measurements means that the system is 
reinitialized in an eigenstate of the measured quantity after each measurement the 
total time needed to take the measurements scales as the square of the number of 
measurements, ignoring the time required for read-out and re-initialisation of the 
qubit . 

We begin by introducing the interleaved sampling (or signal reconstruction) 
technique in section[2] Then in section[3]we show how it can be applied to a signal 
that is generated by a series of projective measurements. Also we describe the qubit 
model that we use to illustrate the benefits of interleaved reconstruction. Section |4] 
presents examples of the reconstruction method and the demonstrates the superior 
accuracy of the interleaved sampling method for estimating the decoherence rate 
of the qubit. In section [5] we conclude. 

2 INTERLEAVED SAMPLING 

The standard approach to the reconstruction of classical signals is to use the 
Nyquist-Shannon sampling theorem to determine the minimum sampling rate. If 
a signal, x(t), has no frequency components above a finite cut-off frequency (S), 
then the signal can be reconstructed exactly by sampling at a rate of fs = 2B 
and then interpolating the signal between the sample points using normalized sine 
functions [7]. For a series of samples x[n\ (where n is an integer, and x[n\ corre- 
sponds to a measurement at time tn = nAt = n/2S), the interpolated function is 
given by 

oo t t 

x{t) = ^ a;[n]sinc (^^) ■ (1) 

n— — oo 

The reconstructed signal is an exact reconstruction of x{t) provided that x{t) is 
limited to the band [0,i?]. In the frequency domain the sampling process causes 
the reconstructed signal to be copied into frequency bands at integer multiples of 
fs ; because the signal is band-limited the original spectrum does not overlap with 
these copies. 

In practice, real signals are not limited to a finite band of frequencies. The 
frequency transform of a causal signal can be zero at discrete values but not 
over any finite band [8]. Thus the signal will not be completely confined to the 
required frequency band. When the signal is not band-limited, significant portions 
of the spectrum outside the relevant frequency band will be copied into the band 
of interest during the sampling process. This is referred to as aliasing. However, 
as long as the signal components outside the band are negligible, the distortion 
introduced by aliasing will be small. 

The minimum sampling rate of fs = 2iJ and the reconstruction formula, 
Eq. ([1]), are applicable to 'baseband' signals (i.e. ones that contain low frequency 
components near DC and up to some maximum frequency B). For signals that 
do not contain significant components at low frequencies, it is possible to reduce 
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the sampling rate by generating two lower-rate measurement series and then use 
a more complex reconstruction formula, due to Kohlenberg [9], which interleaves 
the two series of measurements. If the signal lies between a lower frequency and 
an upper frequency such that /[/ = /l + B, where the signal bandwidth is B, the 
two sample rates for the interleaved sampling need only be fs = B. That is, the 
total sample rate is 2B rather than the fs = 2{fL + B) that would be required 
if the standard interpolation/reconstruction formula were to be used. If ^ B, 
the saving in measurements is considerable. 

The two sampling processes are taken at intervals At = \ / B and separated by 
a time k. The separation between samples, k, can take any value, < fc < At, as 
long as the interpolation function is well-behaved. The generalized interpolation 
formula is given by [10] . 

x{t)= [xo[n]S{t-tr,)+Xk[n]S{-t + t^^k)] (2) 

where the two series a;o[n] and a;fc[n] are measurements at tn = nAt and tn.fe = 
{nAt -\- fc) respectively, and the interpolation function S{t) is 

S{t)=SQ{t) + Si{t) (3) 

where 

_ cos[27r(rg - fL)t - rirBk] - cos[2TTf Lt - rnBk] 

2TTBt8in{rTTBk) ' ^' 

cos[27r(/L + B)t -{r+ l)7rBfc] - cos[27r(rS - fL)t - {r + l)7rBfc] 
~ 27rBtsin[(r-hl)7rBfc] ' 

(5) 

where r is the smallest integer larger than The interpolation function has 
the values: S'(O) = 1 and S{nAt) = S{nAt + k) = 0. For the special case where 
the lower frequency is zero and k = At/2 = 1/25, the interpolation formula 
reduces to 

We now consider how this interleaved sampling method could be used to char- 
acterize the Hamiltonian and decoherence properties of a qubit. 



3 QUBIT MODEL AND MEASUREMENT 

There are several key differences between the reconstruction of a classical signal 
and the reconstruction of quantum coherent oscillations. The primary difference 
is the unavoidable backaction in quantum measurements. Quantum systems can- 
not be continuously measured without disturbing the dynamics one is trying to 
probe. Continuous measurements are possible, however continuous quantum pa- 
rameter and state estimation can be computationally expensive, see Refs. [llj and 
[12| . It is for this reason that most treatments are restricted to periods of evo- 
lution punctuated by projective measurements of an observable az for example. 
Such measurements only ever provide a +1 or —1 measurement result, with the 
probability of the result ±1 determined by the quantum state. The results ±1 
correspond to the |+) = (1,0)^ and |— ) = (0,1)^ eigenstate of the measured 
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quantity. The sample values that we need to reconstruct the signal are found by 
taking a series of measurements - which is termed the record (one realization of 
the random process) - and averaging the record over multiple realizations. There 
are TV measurements at each sample point in the record, giving a measurement 
error (projection noise) which scales as \/y/N . In the limit A'' — > oo, the projection 
noise vanishes. After each measurement the quantum system has been projected 
into one or other of the eigenstates - effectively reinitializing the system - and the 
evolution/oscillations must begin again. This means that the total time taken for 
a series of M sample points is at least 

Tmin = N{At + 2At... + MAt) = \NM{M + I) At (6) 

so the total time taken to generate M sample points scales approximately as the 
square of M , and any reduction in the number of points has a large effect on 
the minimum time taken to characterize the qubit. However, the time interval 
At is longer for the interleaved sampling method because the interleaved sample 
rate is lower than the Nyquist-Shannon sampling rate, so the saving in the total 
measurement time is approximately equal to the ratio of the sample frequencies. 
It is worth noting that there is also a saving in the time required to generate 
the measurement record by not re-preparing the initial (+1) eigenstate after each 
measurement if the —1 result is obtained. This can be achieved in software for 
both the sine and interleaved reconstruction methods. 

The state of the qubit is represented by the density matrix p and its evolution 
can be found by solving the master equation. For a Hermitian Lindblad operator 
c = ^/K,y = c\ where k is the strength of the environmental interaction, this is 
given by [13] 

p = -i H,p 



^[y\y,p\]- (7) 



The typical decoherence timescale is t = 1/k. This equation assumes that the 
environmental noise is uncorrelated (Markovian) and that the coupling to the 
environment is relatively weak. 

We can rewrite the Master equation in the Bloch vector representation [13] . 
with 

= ^Y[a^p] i e {x, y, z} (8) 
where ai is a Pauli matix and the density operator is 

P = ^{I + rxO-x + Tyay + r,az). (9) 

With a Hamiltonian H = ujaxfi (where oj = 27r/ and / is the characteristic qubit 
oscillation frequency) and ij = a^, equation ([7]) provides three coupled equations 

fy = —UJrz — ^KTy, (10) 

rz = i^ry. 

which can be solved analytically. For an initial condition rz = 1 and rx = 0, 

rxit)=0 (11) 
ry(t) = --e-^'"*sin(Ait) (12) 

rz{t) = e"^**^ [cos(t^) + 2k sin(t^)/^] (13) 
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where /i = yj uj^ — Ak'^. Restricting ourselves to a measurement of the z component 
of spin, for simplicity, the probability of getting a +1 measurement result at a time 
t is simply P+{t) = {rz + l)/2 and P-{t) = 1 — P+{t) for the —1 result. We can use 
the analytic expressions to find the probabilities for the result of each measurement 
as a function of time and generate a series of A'' measurements to simulate a set 
of experiments. By generating a set of A'' results at each of the sample points, 
we can average the results to find an estimate the signal level/oscillation at the 
appropriate times and then use these values in the interpolation/reconstruction 
formula ([2]) as if it were a noisy classical signal. 

The Hamiltonian and measurement were selected as an example because they 
represent a common solid state device: a charge qubit quantum dot with a tun- 
neling interaction and a charge measurement [Tl[2l[3l[T3] . 

4 RESULTS 

Figure [1] shows an example of the reconstruction of coherent oscillations, that 
is Eq. (|13p . using both the sine [see Fig. [1] (a)] and interleaving methods [see 
Fig. [1] (b)]. The value of k selected is at the upper end of the range for which 
the technique can be used (in practice k < 0.3). The lower the value of k, the 
narrower the spectral profile and the better the approximation that the signal is 
confined to the selected frequency band. Ideally, the signal outside the frequency 
band should be as low as possible, but there is a trade-off between widening the 
band and reducing the number of measurements used in the interleaving. If the 
decoherence is sufficiently weak for the frequency components outside the band 
[fL,fu] to be small, then the interleaving technique is remarkably robust. From 
figure [T] it is apparent that the sine reconstruction is accurate and appears to 
be robust to projection noise. This must be evaluated with regard to our aims 
of time-efficiently characterizing the fundamental frequency and the decoherence 
parameter. With this in mind consider the following: the sine method is using 4200 
samples compared to 1800 samples for the interleaving method. In this example, 
the number of samples for the interleaved method is approximately 0.43 times the 
number of samples required by the standard sine method. 

Figure 2 shows the Fourier amplitude of the frequency components in and 
around a frequency band fj^ = 0.8 to /[/ = //,+ _B = 1.2 for a central frequency 
of / = 1, K = 0.02 (t » 50), N = 100, and M = 160 (over 200 osciUation cycles). 
The number of measurements required for the interleaved method is three times 
smaller than that required from Nyquist-Shannon sampling, and the minimum 
time required to take the measurements is reduced by a similar factor (x3.03). 
The central peak in the frequency band is clearly visible and the figure also shows 
a Lorentzian resonance curve that has been fitted to the data in the frequency 
band. The main frequency components are present in the band [Jl, fu\ - Figure 2 
also depicts the (noisy) tails of the distribution within the band; there are lim- 
ited components outside this band. There are aliased peaks at higher frequencies, 
around twice the central frequency and above, but the amplitude of these higher 
frequencies is much lower than the main Fourier peak shown in Figure 2. 

Both the central frequency and the width of the resonance (e.g. the full peak 
half maximum) can be found by fitting a resonance curve to the Fourier transform 
of the reconstructed signal. This provides an estimate of the central frequency (and 
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Fig. 1 (Color online) Signal reconstruction with k = 0.1, / = 1, = 0.8, B = 0.4003. 
The dashed lines are a plot of Eq. JTSj. (a) the solid blue line is the sine reconstruction, i.e. 
Eq. 1^ of coherent oscillations from the noisy signal from = 100 projective measurements 
per point (M = 42, two samples not shown). The solid circles represent the averaged record, 
which is plotted at times l/fs (b) the squares and circles are the two interleaved (averaged) 
records with N = 100 (M = 18, one sample not shown). The solid blue line is the interleaved 
reconstruction, ic. Eq. 




Fig. 2 (Color online) The frequency distribution for a reconstructed signal using the inter- 
leaved measurements (blue-solid line), with fitted resonance curve (red-dot-dashed line) and 
the sampling frequency band {B = 0.4, black-dashed lines). 
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Fig. 3 (Color online) Estimates for the decoherence time from the fitted resonance (r = 

for different number of projective measurements per sample point: f = 1, k = 0.02, B = 0.4 

and 200 oscillation cycles, mean values and errors. The error bars are iter. 
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Fig. 4 (Color online) Convergence of the estimated decoherence time (r = 1/k) as a function 
of the number of samples M for the sine method (blue-dash) and the interleaved method 
(black-solid). All points use N = 100 and / = 1, k = 0.05, B = 0.4 . The error bars are ztcr. 
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hence the Hamiltonian) and the decoherence rate, which is related to the width 
of the peak. The accuracy of the frequency estimate is dependent on the spacing 
of the Fourier coefficients and can be improved slightly by interpolating between 
Fourier coefhcients [Tl] . For the cases considered here the frequency errors are less 
than about 1% and the main errors come in the estimation of the decoherence 
rate. The resonance curve corresponding to the coherent oscillations given in (|13p 
is fitted to the signal peak in Figure 2 using a simple enumerative search through 
a range of parameter values (for frequency, k or t and the peak amplitude). 

The average decoherence time (as given by t = 1/k) is plotted (with error 
bars) in Figure 3 against the number of measurements (A'^) for each sample point. 
The mean value contains a statistical bias which underestimates the decoherence 
time (overestimating the environmental coupling) because the effect of aliasing is to 
increase the size of the Fourier coefficients in the tails of the resonance peak, which 
broadens the measured peak. If more measurements are taken for each sample 
point, the representation of the coherent oscillations will be better, the level of 
the noise in the frequency domain will be reduced and the frequency components 
outside the band of interest will be smaller, thereby reducing the effect of aliasing. 
The minimum time required to take the measurements to characterize the qubit 
is only linear in the number of measurements at each sample point A'', and Figure 
3 shows that TV > 50 gives a reasonable accuracy for the decoherence time. 

Figure 4 demonstrates the superiority of the interleaving method over the sine 
construction. For M > 150 the interleaved estimate of t provides a good estimate 
(t = w 50) while the sine method requires significantly more samples {M ^ 
300) to achieve similar results. Similar results were obtained for a large range of 
values for N, M and k with / = 1 fixed. 



5 CONCLUSIONS 

In this paper we have demonstrated that a sample (and time) efficient strategy 
exists for using a series of projective measurements to estimate the frequency of 
the coherent oscillations of a qubit subject to a decohering environment. This 
strategy uses an interleaved sampling technique from classical signal processing 
for the reconstruction of band-limited signals. If the qubit resonance is confined 
to a narrow frequency band, this technique could offer significant savings in the 
numbers of measurements that are required to characterize the system and in 
terms of the length of time required to take the measurements. We have used a 
simulated example to generate sets of sample points that demonstrate that the 
technique is robust even when the coherent oscillations are not strictly confined 
to the relevant frequency band and when there are errors in the sample points 
themselves. 
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